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Abstract 

The accessory beta subunit (Ca^P) of calcium channels first appear in the same genome as Ca^l L-type calcium channels in 
single-celled coanoflagellates. The complexity of this relationship expanded in vertebrates to include four different possible 
CavP subunits (P,, fi2, Ps, P4) which associate with four Ca„l channel isoforms (Cavl.l to Cav1.4) and three Ca^l channel 
isoforms {Cav2.1 to Cav2.3). Here we assess the fundamentally-shared features of the CavP subunit in an invertebrate model 
(pond snail Lymnaea stagnalis) that bears only three homologous genes: {LCa^^, LCb^,!, and LCa^P). Invertebrate Ca^P 
subunits (in flatworms, snails, squid and honeybees) slow the inactivation kinetics of Ca^l channels, and they do so with 
variable N-termini and lacking the canonical palmitoylation residues of the vertebrate P2a subunit. Alternative splicing of 
exon 7 of the HOOK domain is a primary determinant of a slow inactivation kinetics imparted by the invertebrate LCBvP 
subunit. LCa„p will also slow the inactivation kinetics of LCav3 T-type channels, but this is likely not physiologically relevant 
in vivo. Variable N-termini have little influence on the voltage-dependent inactivation kinetics of differing invertebrate CavP 
subunits, but the expression pattern of N-terminal splice isoforms appears to be highly tissue specific. Molluscan LCa^P 
subunits have an N-terminal "A" isoform (coded by exons: la and lb) that structurally resembles the muscle specific variant 
of vertebrate pi a subunit, and has a broad mRNA expression profile in brain, heart, muscle and glands. A more variable "B" 
N-terminus (exon 2) in the exon position of mammalian P3 and has a more brain-centric mRNA expression pattern. Lastly, 
we suggest that the facilitation of closed-state inactivation (e.g. observed in Cav2.2 and Ca^Ps subunit combinations) is a 
specialization in vertebrates, because neither snail subunit (LCav2 nor LCavP) appears to be compatible with this observed 
property. 



Citation: Dawson TF, Boone AN, Senatore A, Piticaru J, Thiyagalingam S, et al. (2014) Gene Splicing of an Invertebrate Beta Subunit (LCavP) in the N-Terminal and 
HOOK Domains and Its Regulation of LCavl and LCav2 Calcium Channels. PLoS ONE 9(4): e92941. doi:10.1371/journal.pone.0092941 

Editor: Alexander G. Obukhov, Indiana University School of Medicine, United States of America 

Received April 17, 2013; Accepted February 27, 2014; Published April 1, 2014 

Copyright: © 2014 Dawson et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits 
unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. 

Funding: This work was funded through a Heart and Stroke Foundation of Canada Grant-ln-Aid and NSERC Discovery Grant to J.D.S. and an NSERC Canada 
Graduate Scholarship and Ontario Graduate Scholarship award to A.S. The funders had no role in study design, data collection and analysis, decision to publish, or 
preparation of the manuscript. 

Competing Interests: J. David Spafford is a PLOS ONE Editorial Board member. He also confirms that his Editorial Board membership with PLOS ONE does not 
alter his adherence to all the PLOS ONE policies on sharing data and materials, as detailed online in the PLOS guide for authors. 

* E-mail: spafford@uwaterloo.ca 



Introduction 

Unique ancillary beta subunits are ideirtifiable in the proteo- 
somal complex with different voltage-gated ion channels including 
K [1-3], Na^ [4~6] and Ca [7-9] channels. These accessory 
subunits are known to promote the membrane expression and 
trafficking of ion channel complexes, as well as to modify the 
biophysical features of voltage-gated ion channels [7-9] . Sodium 
channels (Na^,l.x) and calcium channels (Ca^,l.x, Ca^,2.x, Ca^,3.x) 
bear a common structural template of four repeat domains of six 
transmembrane helices each, but they have distinct ancillary 
subunits which are known to regulate them [7-9] . Sodium channel 
beta subunits (Na^P) evolved separately in different animal groups 
such as insects, snails and vertebrates. TipE/Teh are insect Na„P 
subunits and have a likeness to (Slo/BK) KvP subunits with EGF- 
like domains [4], while gastropod snails Na„P subunits is a CUB 
domain containing protein family [10], while vertebrate Na„P 
subunits are CaM-like with a V-set Ig extracellular loop [5]. 



Calcium channel beta subunits (Ca^P) have common homologs in 
animal groups rooted in genomes of single cell organisms like 
coanoflagellates, which also possess an L-type calcium channel 
(Ca^l) homolog [11]. From a likely primordial template of one 
Ca^l and Ca^P subunit in coanoflagellates [12], emerged the 
complexity in numbers in vertebrates, where there are four 
different Ca^,P subunits, (Pi, P2, P3, P4), which are expected to 
regulate seven pore-forming oti subunits from the two high 
voltage-activated classes of calcium channels, Ca^Ll-Cavl-^ (L- 
type) and Cav2.1-Cav2.3 (non-L-type) [13]. There are also low 
voltage-activated calcium channels, Cav3.1-Cav3.3 (T-type), but 
there is not strong evidence for their regulation by Ca^P subunits 
[7]. 

The pond snail Lymnaea stagnalis, possess a simple set of calcium 
channels: a single L-type channel gene (LCa^,l) [14,15], a single, 
synaptic, non-L-type channel (LCa^2) [16,17], a T-type channel 
(LCa„3) [18,19] and a single beta subunit, LCa^,p [20]. We use this 
simple model to address what are the core conserved features for P 
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subunit regulation of calcium channels, and to determine what 
features are adaptive and unique to vertebrates. 

We examine the common alternative splicing patterns in the N- 
terminus and HOOK domains between snail LCa^.p and 
vertebrates CavP subunits. mRNA expression suggests a tissue 
specificity to the expression of N-terminal isoforms, with LCa^,PA 
(eg. structurally resembling vertebrate pi isoform in sequence) 
possessing a more generalized expression, while the LCavPB (in 
exon position of the shorter vertebrate p,3 isoforms) has a more 
discrete, brain specific expression pattern. Molluscs possess a 
common site of alternative splicing in the HOOK domain, which 
generates a variably sized exon 7 by use of alternative acceptor 
site, and is in the same location as mutually-exclusive spficed 
isoforms (exon 7a,7b, or 7c) [8,9] in vertebrate Ca^P subunits. 
Ca^P subunits in invertebrates (eg. schistosomes [21], snails [20], 
squid [22] and honeybees [23]) commonly slow the inactivation 
kinetics of Cav2 channels in a manner like vertebrate P2a subunit 
[24-26], but without the canonical N-terminal palmitoylation site 
that primarily influences its slow inactivation. We show that the 
HOOK domain (- isoform) more than the N-terminus (A isoform) 
of snail LCa^P contributes to the slowing of inactivation kinetics of 
Ca^2 channels. The snail Ca„P subunits can also slow the 
inactivation of snail LCav3 T-typ(- channels, but this may be an 
observable phenomenon in vitro, because there is littie evidence for 
CavP subunits associating with T-type channels in native cells [7]. 
And lasdy, we examine that the capacity that CavP subunits to 
promote an increase in the inactivation rate of 0( subunits while the 
channels are largely in the closed state [27,28]. This property 
known as "closed-state inactivation" is not observable for the snail 
LCavP or synaptic LCav2 channel, which suggests that this 
property is more a specialization for specific vertebrate subunit 
combinations such as CbvPs and Cav2.2 [27,28]. 

Results and Discussion 

1. Conserved SH3-GK core of CavP subunits 

P subunits of calcium channels (Ca^P) contain four major 
representatives in vertebrates (Pi to P4), represented by a single 
gene in invertebrates (Figure lA). CavP subunits are unique 
members within the MAGUK (membrane associated guanylate 
kinase) family of proteins [1 1,29,30], which share a structural core 
consisting of an SH3 domain and a guanylate kinase (GK) domain, 
separate by a variable HOOK region (Figure. IB, IC). 
MAGUKs, are noted for their role as scaffolding and cytoskele- 
tal-organizing proteins, such as post-synaptic density protein 
(PSD-95), where the SH3 and GK domains form potential 
intramolecular or intermolecular interactions, with the same or 
difiering MAGUKs, with neither the SH3 nor GK domains have 
active sites for standard polyproline interactions or nucleotide 
kinase activity respectively [11,29,30]. P subunits of calcium 
channels lack PDZ domains of standard MAGUKs, but does have 
a HOOK domain which splits the SH3 domain [31—33] 
(Figure IB, IC). 

The alpha helix at the end of hydrophobic segment 6 in 
Domain I of Cavl and Cav2 ai subunits, is largely expected to 
extend to an alpha 1 -interaction domain (AID) of conserved amino 
acids (blackened amino acids, Figure ID), which in a crystal 
structure embeds into a deep hydrophilic groove of the CavP 
subunit formed by the guanylate kinase domain [31-33] (black- 
ened amino acids. Figure IC). The SH3 and GK domains of 
CavP subunits (Figure IB), and most notably the key-in-keyhole 
interaction residues, between the AID sequences of calcium 
channel oel subunits and a conserved binding pocket formed by 
the CavP subunit GK domain, respectively, are largely conserved 



down to the simplest known organisms that have calcium 
channels, the unicellular eukaryotes choanoflagellates 
(Figure IC, ID). The most basal unicellular organisms (proto- 
zoan) or multicellular animals (sponge, poriferans) have a single 
Csvl L-type calcium channel homolog (Figure lA). These Cavl 
homologs, in the simplest organisms, are structural predecessors to 
the high voltage-activated, dihydropyridine-sensitive LCavl ho- 
molog expressed in the pond snail, Lymnaea stagmlis [14,15]. A 
second and third class of calcium channels likely derived from L- 
type calcium channels (in ancestral relatives of placozoans, 
cnidarians) include the non-L-type, Cav2 channel class, noted for 
their unique roles in mediating synaptic transmission [16,17], and 
the low voltage-activated T-type, Cav3 channel class [18,19]. AU 
Cavl and Cav2 channels from simple representatives (placozoan, 
cnidarians) have a hallmark AID sequence [34] for associating 
with CavP subunits, which enables a promiscuity and interchange- 
ability of CavP subunits interacting with difiFering Cavl and Cav2 
ai subunit classes. 

2. Variable regions are associated with common 
alternative splicing 

Variable regions outside of the conserved SH3 and GK 
domains of CcvP subunits provide specificity and unique 
modulation of of different Cavl and Cav2 channel types 
[8,9,13]. CavP subunits dramatically vary in the N-terminus 
(VI), HOOK region (V2) and C-terminus (V3) which is illustrated 
in a running similarity score (Figure IB) of aligned sequences of 
invertebrate and mammalian CavP subunits (Figure SI). Alter- 
native splicing of CavP subunits coincides with these variable 
regions, illustrated in the exon-intron structure of snail LCavP and 
the four human CavP isoforms (shown in Figure 2A). 

The major pattern of alternative splicing in the four vertebrate 
CavP subunits is represented in the moUuscan LCavP subunit (as 
shown in Figure 2A), suggesting that the gene splicing patterns 
evolved in the common ancestor of moUusks and vertebrates, 
before the duplication that lead to the generation of the our 
vertebrate isoforms. 

3. Conservation of N-terminal mutually-exclusive splicing 
in CavP subunits (exons la/lb and exon 2) and their 
consistent large intron sizes 

The two most upstream CavP subunit exons form one N- 
terminal splice variant, dubbed variant CavPA (Figure 3A), 

consisting of exon la with highly variable sequence, jointed to 
exon 2b which contains a common conserved amino acid clusters 
(KxSDSG...FlRQ} shared between moUuscan CavPA and verte- 
brate spliced isoforms pi, P2c, P2d, P4b. The mutually exclusive 
N-terminal spUce variant CavPs is formed by exon 2 (Figure 3B) 
and is downstream of exon la and exon lb in the genomic 
structure of CBvP subunits (Figure 2A). CavPs generates short N- 
terminal vertebrate isoforms P2b, P2e, P3, P4a including the 
canonical P2a which has doublet cysteines that are palmitoylated 
and confer slow inactivation to vertebrate Cavl and Cav2 calcium 
channels [24-26], and to moUuscan LCav2 channels [20]. 
Honeybee, as a representative insect and arthropod, has two N- 
terminal spliced isoforms [23], but does not possess a conserved 
exon IB with the (KxSDGS...FIRQ) cluster found in moUuscan 
[22] or vertebrate [8,9] CavPs isoforms (Figure 3). 

A general feature of N-terminal exons la, lb and 2 is that the 
intervening non-coding sequences (introns) that span the region 
are so large that they compose the majority of the genomic region 
spanning most CavP subunits even in a simple cnidarian species, 
Msmatostella (Figure 2B). The intron sizes are largest in vertebrates 
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Figure 1. Conservation of calcium channel beta (CavP) subunits from single-celled organisms to humans. (A) Gene tree derived from 
aligned sequences illustrating the relationship of the four vertebrate Ca„p subunits ((31, )32, [33, p4) with the single Cavp representatives in non- 
vertebrates, including single-celled organisms. (B) Running average of similarity (window = 50 amino acids) of aligned invertebrate and human Ca„p 
subunits sequences (in Figure SI). Cavp subunits have highly conserved SH3 and guanylate kinase (GK) domains, with variability in the N-terminus 
(VI ), HOOK domain (V2) and C-terminus (V3). Alternative splicing shared in CavP subunits in Exons 1 a/1 b and exon 2 (N-terminus), and exon 7 (HOOK 
domain). (C) IVlultiple alignment of Cavp subunits from single cell organisms to humans illustrating highly conserved SH3 and GK domains, conserved 
secondary structures (a helices and CavP sheets), and calcium channel (AID) binding residues (blackened residues) reported in crystal structures of 
CavP subunits [22-24]. Exon boundaries are indicated in red. (D) Alignment of cytoplasmic region post trans-membrane segment 6 in Domain I of 
Cavl and Cav2 calcium channels which forms an expected a helix. Regions include a proximal linker (PL) and the Alpha, -Interaction Domain (AID) 
associated with binding CavP subunits (blackened residues). 
doi:1 0.1 371 /journal.pone.0092941 .gOOl 



isoforms (such as fish and mammalian subunits [35,36]) which also 
have the largest number of exon variants that code for exons 1 a, 
and exon 2 in the N-terminus, which are 6 and 4 exons 



respectively, in vertebrate CavP2 and CavP4 subunits (Figure 2A, 
2B). The genomic region spanning the N-termini for hCavP.i and 
hCavP2-is 217,000 bp and 236,000 bp. Gargantuan introns such 
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Figure 2. Conserved exon-intron organization, alternative splicing, and N-terminal intron sizes in the genomic sequence spanning 
calcium channel beta (Cavp) subunits. (A) Alignment of the 15 exons of Cav|3 subunits comparing snail and human Ca^P subunit splicing. Ca^P 
subunits have mutually exclusive splicing of N-terminal exon la/lb or exon 2 isoforms. Exon 7 in the HOOK domain is subject to mutually-exclusive 
exon splicing (exon 7a or exon 7b or exon 7c) in vertebrates or splicing in mollusks generated by alternative acceptor sites (exon 7a+, exon 7a-). 
Molluscan and vertebrate have truncated forms of CavP subunits lacking the GK domain and C-terminus as a result of skipping of exon 7. (B) Most of 
the intron sizes of CavP subunits span the N-terminal exons, and the size of the total intron sequence in the N-terminus increases with the number of 
exons in the N-termini. 
doi:l 0.1 371/journal.pone.0092941 .g002 



as in the N-termini of hC£vP4 and hCavPs are highly unique, and 

represent ~0.1% of all introns in the human genome [37]. The 
evolutionarily conserved, larger sized introns are expected to 
contain DNA elements which regulate the timing and tissue 
specificity of the expression of GOvP subunit isoforms. 

4. mRNA expression confirms a tissue specific expression 
pattern for N-terminal A and B isoforms 

mRNA expression patterns suggest that the N-terminus of 
LCa„P, specifically, exons la/lb (dubbed "LCavPA") and exon 2 
(LCa„Pj3") are highly regulated in their tissue expression 
(Figure 4A), while the HOOK domain splice variants are not 
(Figure 4B). LCa^PA is a widely expressed isoform in most tissues 
(Figure 4A, top), and resembles the expression profile of the 
HOOK domain splice isoforms which lack a tissue specific 
distribution (Figure 4B). In contrast, LCa^PB is more discretely 
expressed than LCa„PA, as a mostly brain specific isoform lJ2°/o), 
with residual expression levels in the heart (16%) (Figure 4A, 
bottom). Comparing the relative mRNA levels of the Ca„P 
(Figure 4A) with the corresponding mRNA expression of LCayl 
and LCav2 channels in different snail tissues (Figure 4C) 
resembles similarities in expression patterns to mammalian gene 
spliced isoforms [8,9]. LCa^PA (Figure 4A) is the primary Ca^P 
subunit isoform expressed in the equivalent tissue resembling 
skeletal muscle in snails (buccal mass and foot), where LCayl 
channel expression dominates (Figure 4C). In contrast, LCa^PB 
and LCav2 channel are almost undetectable in the equivalent 
tissue of skeletal muscle in snails (Figure 4A, Figure 4C). The 
combination of mostly LCa,,l and LCa^PA in snail muscle is 
consistent with the exclusive pairing of Ca^Ll and CavPla in 
vertebrate skeletal muscle [38-40], where vertebrate Ca^pi bears 
only N-terminal exons la/lb like snaU LCavPA (Figure 2A). It is 
also noteworthy that mammalian P3 lacks N-terminal exons la/ lb 
and is composed of exon 2 only (Figure 2A) like the snail LCa^Pg 
isoform. Mammalian Ca^Ps is mostiy a brain specific subunit that 
pairs more often with synaptic Cav2.2 channels [41,42], and in this 
regard, resembles the sphce variant of snaU LCavPs which is also is 
more brain specific variant (Figure 4A) and composed of exon 2. 
These are examples of likely common pairings of calcium channel 
and beta subunit isoform between invertebrates and vertebrates, 
but outside of this, there appears to be a lot of promiscuity and 
overlap between Ca^l and Cav2 channels and their association 
with the differing CavP subunit isoforms. It is noteworthy that "A" 
isoforms are found in examples of all four vertebrate CavP 
subunits, and only CavP^ is lacking a "B" isoform in vertebrates 
(Figure 4C). 

Changes in mRNA levels of CevP subunits from juvenile to 
adult snails, match the direction of change of the corresponding 
mRNA of calcium channels (Figure 4D). Heart and muscle 
mRNA levels fall from juvenile to adult animals for LCavl and 
LCav2 channels and LCavP subunits, while there IS a rise in LCayl 
and LCav2 channels and LCavP subunits corresponding with the 
maturation of sexual organs from juvenile to adult animals, such as 
in prostate and albumen gland (Figure 4D). 



5. Common HOOK domain splicing of exon 7 in 
molluscan and vertebrate Ca^p subunits 

The HOOK domain is a second area of variability in CavP 
subunits and subject to alternative splicing. The HOOK domain is 
unstructured and not resolvable in the crystal structure of CavP 
subunits [31-33]. The HOOK domain splits the SH3 domain, 
separating the 5* CavP strand from the rest of the SH3 domain 
(Figure IB, Figure IC and Figure SI). The variability in the 
HOOK domain is largely contained in exons 6 and 7, where exon 
7 is subject to alternative splicing that is shared in molluscan and 
mammalian CavP subunits (Figure 2A). 

Mutually-exclusive splicing of exon 7 in mammals leads to 
either a skipping of exon 7 altogether, or includes a choice of one 
of three exon 7 variants including, exon 7a (a long form: e.g. 
AIDID...PFFKK), exon 7b (a short form: AKQKQKQ) and 
more rarely exon 7c (Figure 3B) [8,9]. Snails (and also squid [22]) 
can also exclude exon 7 altogether, via use of alternative acceptor 
sites (Figure 2A) or generating exon 7 containing seven or five 
extra amino acids (exon 7-I-), respectively, or versions lacking these 
extra amino acids (exon 7—) (Figure 3B, Figure S2). Inverte- 
brate exon 7 resembles the exon 7a (long form) of vertebrate CavP 
subunits, with a conserved polybasic region at its 3' end 
(KRxPFFKK) (Figure 3B) [43]. The skipping of exon 7 is found 
in mammalian CavP subunits [44—46] and shared in molluscan 
CavP subunits and generates a frame shift and an immediately 
truncated CavP subunit that lacks the GK domain and C-terminus 
(Fig^e 2A). The frame shift occurs because exon 6 ends within a 
codon after the first nucleotide (phase 1), whereas the intron 
preceding exon 8 is located between codons (phase 0) (Figure S2). 
One role discovered for the truncated CavP subunit in vertebrates 
is that it transports into the nucleus and serve as a transcription 
factor [47-50]. Conservation of splicing that skips exon 7 suggests 
that truncated CavP subunits in invertebrates may also have 
similar roles outside of their association with Cavl and Cav2 
calcium channels. 

6. HOOK domain spliced variants lack a tissue specific 
expression pattern 

Snail LCavP+ or LCavP — containing or not containing the 
extension to exon 7, respectively, show no obvious tissue 
preference in their mRNA expression. Notably LCavP — (lacking 
the extension to exon 7) is more common, being 74%— 84% of the 
total mRNA transcript in all snail tissues (Figure 4B). The lack of 
tissue regulation for splicing in the HOOK domain contrasts with 
the more highly regulated N-terminal alternatively-spliced exons. 
There is also consistentiy smaller and more normal sizes of introns 
spanning exon 7, which contrasts with the very large to 
gargantuan introns spanning the N-terminal exons (Figure 2B). 

7. Invertebrate CavP subunits slow the inactivation 
kinetics of Cav2 channels 

We have co-expressed full-length isoforms of LCavl, LCav2 and 
LCav3 channels with mammalian o.2?>\, and the four combinations 
of LCavP isoforms in HEK-293T cells to assess the consequences 
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Figure 3. Multiple sequence alignments illustrate conserved splicing in (A) the N-terminus and (B) HOOK domain of invertebrates 
and vertebrate CavP subunits. Alternate N-terminal isoforms composed of exons la-lb (CavPA) or exon 2 (Ca„PB), and HOOK domain splicing 
includes optional short addendum to exon 7 (Ca„p— /Ca„p+) in invertebrates or mutually exclusive splicing, exon 7a, 7b, 7c. Note that Ca^p from 
snails, squid, schistosomes and bees and vertebrate CavP2a have slow inactivation kinetics and (B) possess HOOK domains with a long form of exon 7 
(A form) with a common polybasic region at its 3' end. 
doi:1 0.1 371 /journal.pone.0092941 .g003 



of LCa^,p splice isoforms on the expression of the calcium channels 
recorded using whole-cell patch clamp electrophysiology. Combi- 
nations that we assessed were LCa^P with exon la/ lb (A form) or 
exon 2 (B form), and LCa^,p channels with or without the 
optional seven amino acids in exon 7 (i.e. -I- form or — form, 
respectively): (LCavPA+, LCa^PB-H, LCa„PA— , LCavPs") com- 
pared to the absence of co-expressed Ca„P subunit. We substituted 



native external calcium ions for barium ions to evaluate the 
consequences of calcium-independent effects on biophysical 
properties recorded in whole cell voltage clamp, avoiding the 
dramatic calcium dependent inactivation, characteristic of LCa^ 1 
channels, observed when external calcium in the charge carrier. 
Reported effects of Ca^P subunits have focused on the dramatic 
changes to the voltage-dependent properties [7-9], so we have 
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Figure 4. Quantitative RT-PCR results show that N-terminal splice isoforms (LCbvPa arid LCavPe), but not HOOK domain splice 
isoforms (LCavp— and LCavp+) of snail Cavp subunits have tissue specific mRNA expression patterns. mRNA levels are illustrated as fold 
change relative to HPRT mRNA levels. (A) More generalized pattern of splicing of LCavpA containing exons la/lb, than LCavPs containing exon 2. 
Exon 2 containing isoform is mostly expressed in the brain, and residual levels in the heart. (B) Exon 7 splicing generates seven extra amino acids 
(exon 7a— vs exon 7a+) and appears to have no tissue selectivity pattern of expression. (B, inset) Percent of LCavp— vs LCavp+, illustrating that 74%- 
84% of all transcripts lack the extra amino acids in exon 7. (C) LCa„1 L-type channel has a more generalized mRNA expression pattern as LCavPA while 
more nervous system specific LCa„2 has a more discrete expression pattern as LCavPs isoforms. (D) Rises and falls in the relative fold changes in 
mRNA levels from juvenile to adult animals are correlated between LCavp, LCa„1 and LCa„2 channel subunits. 
doi:1 0.1 371/journal.pone.0092941 .g004 



limited our study to an examination of barium currents. AH 
electrophysiology results +/— s.e.m. are tabulated in Table 1. 

We found that none of the snail LCa^P isoforms influenced the 
current voltage-relationships such as the threshold voltage at 
which barium currents are beginning to be visible or the voltage of 
peak current for LCavl or LCav2 channels (Figure 5A). This 
absence of a difference is reflected in the activation curves fitted 
with a Boltzmann equation (Figure 5B). Similarly there were also 
no effects of LCa^.P subunits on steady-state inactivation 
(Figure 5C). This is different from the dramatic effects of 
mammalian Ca^,P subunits on mammalian Ca^l and Ca^2 
channels [7-9]. Mammalian Ca^P subunits impart large shifts in 



the voltage-dependence of activation to more hyperpolarized 
voltages by approximately —10 mV to —20 mV [7,8]. Mamma- 
lian Ca^P subunits, except for Ca^P2a also cause an approximately 
10 mV hyperpolarizing shift in steady-state inactivation curves, 
and accelerate inactivation kinetics to varying degrees [7-9]. 
Ca^Pja has doublet cysteines that are palmitoylated in an exon 2 
type N-terminus, which promote membrane association of Ca„P 
subunits [24—26]. Palmitoylation is considered to anchor and 
immobilize CavP2a) causing calcium channels to be more reluctant 
to inactivate (shifting inactivation curves in the depolarizing 
direction), and is responsible for much of the slowing of its 
inactivation kinetics [24—26]. Snail LCa^P resembles Ca„P2a in 
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Table 1. Summary of electrophysiology parameters for Figures 5-8. 





Figure 5 and 7. Voltage-sensitivities of LCa„1, LCa„2, LCa^B cliannels with and witKiout LCa^P splice isoforms 
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Figure 5 and 7. LCa„1, LCa„2 and LCa^B current densities with and without LCa 
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Figure 6 and 7. Quantification of Inactivation kinetics of LCa„1, LCa„2 and LCav2 currents with and without LCavP splice isoforms 
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Figure 8. Analyses of closed-state inactivation for Cav2.2 and LCav2 with snail and mammalian beta subunit combinations 
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promoting the .slowing of inactivation kinetics of LCav2 channels 
(Figure 6B, 6C), albeit to a lesser degree than the slowing of 
inactivation kinetics of LCav2 induced by Ca^Pja (see figure 3e in 
[20]). But notably, the slowing of inactivation kinetics is a 



consistent feature in other invertebrate Ca^P subunits (such as 
schistosome [21], squid [22] and honeybee [23] Ca^fi subunits), 
and these subunits, including snail, lack N-terminal palymitoyla- 
tion residues in the N-termini of their Ca„P subunits. 
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Figure 5. Snail beta subunit (LCBvP) splice isoforms containing either N-terminal exons 1a/1b or exon 2 (A or B) or insert in exon 7 
(7a— vs 7a-i-), boost the membrane expression of snail LCavl or LCav2 channels but has no effect on their voltage-sensitivities when 
external barium is the charge carrier. (A) Current-voltage relationships (curve-fitted with Ohmic-Boltzmann equation) (B) Activation and (C) 
Steady-state inactivation curves (curve-fitted with Boltzmann equation). (D) Current densities (pA/pF) shown as box plot (box = s.e.m., whisker 
range = minimum/maximum values). 
doi:1 0.1 371/journal.pone.0092941 .g005 
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Figure 6. Snail beta subunit (LCavp) splice isoforms containing 
either N-terminal exons 1a/1b or exon 2 (A or B) or insert in 
exon 7 (7a— vs 7a+), slow the inactivation kinetics of LCav2 
channels. Four representative and normalized peak barium current 
traces shown as mean, s.e.m. for (A) LCa^l and (B) LCav2 channels co- 
expressed with LCavP splice isoforms or no Cav|3. (C) Rate of 
inactivation decay reflected in the fraction of maximal peak current at 
250 ms time point. All LCavP splice isoforms like vertebrate Cav)32a slow 
inactivation kinetics of ICayl channels. The slowing of inactivation 
kinetics is maximized with the exon 7a- configuration. 
doi:1 0.1 371 /journal.pone.0092941 .g006 



8. The HOOK domain more than the N-terminus 
contributes to the slowing of the inactivation kinetics 
imparted by invertebrate Ca^P subunits 

We observe that splicing in the HOOK domain of snail 
LCa„|3— lacking extra residues in exon 7, has a more dramatic 
effect on the slowing of the inactivation kinetics of LCav2 channels 
compared to snail LCavP-l- possessing extra residues in exon 7 
(Figure 6B, 6C). In contrast, the differing N-termini of LCa^P do 
not influence the inactivation kinetics of LCav2 (Figure 6B, 6C). 
We observed no measureable electrophysiological differences 
between the expression of LCb^Pa (exon la/ lb) and LCa„PB 
(exon 2) isoforms, suggesting that the N-terminal splicing is more 
associated with tissue localization than to altering biophysical 
properties of snail Ca^l and Ca^2 channels. Size of the N-terminus 
may also be a relevant parameter for LCa„PA and LCa„PB, since 
these variants have similarly sized N-termini (50 vs. 47 amino acids 
respectively). Differing lengths of N-termini of mammalian Ca^Pib 
[51] and Ca^Pa [52] channels are known to correlate with altering 
rates of inactivation in mammalian Ca^,2 channels. 

Comparing different invertebrate Ca^,p subunits (schistosome 
[21], snail [20], squid [22] and honeybee [23]), the N-termini are 
highly variable in sequence, yet all impart slow inactivation 
kinetics of Ca^Z channels in the absence of palmitoylation residues 
found in vertebrate CavP2a Likely more relevant than the N- 
terminus for imparting slow inactivation kinetics is the HOOK 
domain which in invertebrate Ca^P subunits resembles Ca^P2a 
(longer spliced isoform, Exon 7A) with a conserved polybasic 
residue stretch at the 3' end [43,53]. The polybasic residue stretch 
is lacking in vertebrate Pu„ Ps and P4 subunits, containing the 
shorter spliced isoform (Exon 7B instead of Exon 7A). These other 
channel types (Ca„Pii-„ Ca^Pa and Ca„P4) bear faster inactivation 
kinetics than Ca^Pza [8,9] . The conserved polybasic residue stretch 
in Exon 7A contributes to the slow inactivation kinetic phenotype, 
even in the absence of palmitoylation residues in the N-terminus 
found in vertebrate Ca^Paa channels [43,53]. The simplest 
evolutionary hypothesis is that vertebrate Ca„P2a retained the 
extended Exon 7A with the polybasic region in the HOOK 
domain shared in a common ancestral CavP subunit (resembling 
the invertebrate Ca„P subunit) to complement its very slow 
inactivation phenotype governed largely by its unique N-terminus 
of palmitoylation residues. Every invertebrate Ca„P subunit (eg. 
cnidarian, nematode, schistosome, moUusk, insect) resembles Exon 
7A with the polybasic region in the HOOK domain, which may 
suggest that the shorter HOOK domain of Exon 7B was a 
vertebrate adaptation to promote a faster inactivation phenotype 
for CavPih, Ca^Pa and Ca^Pi subunits. 

9. Invertebrate Cavp subunits do not alter the 
inactivation kinetics of LCa^l L-type channels 

Snail LCa^.l channels possess little voltage-dependent inactiva- 
tion without Ca^P subunit expression, so it is not surprising that we 
did not observe additional slowing of inactivation kinetics 
imparted by LCa^,P on LCa^l channels when barium is the 
external charge carrier (Figure 6A, 6C). Snail LCa„l channels 
have a prominent calcium dependent inactivation [14,54] as 
mammalian Cavl.2 channels [55], which has been attributabed to 
a conserved C-terminal IQ, motif and a N-terminal NSCATE 
motif for binding calmodulin in the presence of external calcium 
[14]. The much faster voltage-dependent inactivation of vertebrate 
Cav2 channels has been attributed to the more helical, and rigid 
proximal linker (PL) between the inactivation gate of the 
transmembrane segment 6 (S6) helix of domain I and the AID 
sequence for Ca^.p subunit binding [56-58] . Both invertebrate and 
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vertebrate Covl channels have a conserved glycine in the proximal 
linker (Figure ID) which is expected to lower the alpha-helicity 
and rigidity of the helix, and weaken the observed influences of the 
Ca^P subunit on the inactivation kinetics of Ca^l channels [57]. 
The conservation of mammalian calcium channel features in the 
most primitive multicellular animals is consistent with the notion 
that a common ancestor before the appearance of nervous systems 
(i.e. resembling an extant placozoan), possessed a Cav2 channel 
homolog with a strongly voltage-dependent regulation by Ca„P 
subunits, and a Cavl channel homolog with a calcium-dependent 
regulation of inactivation kinetics by calmodulin. 

10. Invertebrate CSvP subunits slow the inactivation 

kinetics of Cav3 channels in vitro 

There have been reports that Ca^S T-Type channels can be 
modulated by the co-expression of CavP subunits in vitro [59-61]. 
We report that the snail CavP can promote the slowing of 
inactivation kinetics of LCayS channels in a manner similar to 
hCsL^2 channels (Figure 7E,F), and this modulation occurs 
without any effect on other biophysical properties 
(Figure 7A,B,C) or membrane expression levels (Figure 7E). 
This commonly observed phenomenon could be an artifact of in 
vitro studies, since there isn't strong evidence that Ca^S channels 
are modulated with Ca„P subunits in vivo. Reported P subunit 
interactions are weak with T-type channels [59-61], and require 
the transfer of 1-11 linker sequences from Cavl or Cav2 channels 
to acquire characteristic features of CavP subunit modulation [62]. 
Also, native Cav3-CavP protein complexes have never been 
identified, knockdown of CavP subunits in native cells do not 
alter T-type channel properties [63.64], and snail LCav3 [18,19] 
and mammalian [65] T-t\pe channels arc easily reconstituted and 
resemble the features of native currents in vitro, without having to 
co-express CavP subunits. 

11. A lack of closed state inactivation in the Ca^P subunit 
of invertebrates 

Another feature that has been attributed to CavP is the tendency 
of Cav2 channels to inactivate more than expected during a train 
of action potentials [27,28]. This is a feature dubbed "closed-state 
inactivation", where CavP subunits promote a greater inactivation 
during the repolarization period, when Cav2 channels are 
expected to be more unavailable, in non - open and inactivated 
states [27,28]. Observation of this persistent inactivation can be 
obtained by using a two-pulse protocol, where a prepulse step 
depolarization is applied to generate maximal currents, and then 
the percent of recovery of the same peak currents assessed by a 
second voltage step after a time delay (e.g. 0.5 to 40 ms) [28]. With 
a short time delay between voltage pulses, such as 8 ms, there is a 
"dip" in the time of recovery observed for mammaUan Cav2.2 
channels co-expressed with CavP^, in particular where a greater 
than expected fraction of Cav2.2 channels are observed to be 
unavailable for opening compared to more brief time periods 
(0.5 ms) or longer time periods (40 ms) of inactivation recovery 
between the prepulse and the test pulse (Figure 8A, bottom 
panel). We observe that the closed-state inactivation is a property 
that is most evident for CavPs and less for CavPn,, and not 
observable for CavP2a (Figure 8A, bottom panel), as has been 
reported previously [28]. Snail LCavPA+ is similar to CavP2a in 
lacking closed state inactivation (Figure 8, sample trace 
shown in upper panel). Not only was the snail CavP unable 
to bestow a closed state inactivation on mammalian Cav2.2 
channels, but the snail LCav2 channel was not compatible with 
closed state inactivation either. None of the CavP subunits 



including LCavPA+ or mammalian CavPa or CavP2a could 
generate the closed state inactivation channel behavior on LCav2 
(Figure 8B). 

12. Calcium current density changes due to CavP subunit 
expression 

CavP subunits are required for expression of Cavl and Cav2 
channels and little or no surface expression in vitro and in vivo is 
observed in their absence. Binding of CavP subunits to calcium 
channels in the Alpha 1 Interaction Domain (AID) sequence of the 
hnker is modeled to promote membrane trafficking and expression 
of calcium channels, by facilitating their export from the 
endoplasmic reticulum [66] and preventing ubiquitination and 
proteosomal degradation [67-69]. Previously we were unable to 
measure an increase in current density of expressed snail LCav2 
channels with co-expressed snail LCavP, when LCav2 channels 
contained an N-terminus of Cav2.2 [20]. We show here that 
expression of both LCayl and LCav2 with its native N-terminus 
produced statistically significant increases in current density with 
the expression of all snail LCavP subunit isoforms (boxplots in 
Figure 5D, and see statistics in Table 1). 

13. Conclusions 

CavP subunits found in single-celled coanoflageUates to humans, 
bear highly conserved functional core of SH3 and GK domains for 
associating and modulate high voltage-gated calcium channels. 
Variable regions in the N-terminus (exons 1 a/ 1 b and 2) and 
HOOK domain (exon 7) contain conserved alternative splicing 
that is shared amongst CavP subunits. MoUuscan (squid [22], snail 
[20]) schistosome [21] and insect (bee) [23] CavP subunits all 
substantially slow inactivation kinetics of Cav2 channels despite 
differing spliced N-termini of exon la/ lb or exon 2 isoforms. The 
slowing of inactivation kinetics of invertebrate CavP subunits 
occurs in the absence of the canonical paknitoylated residues 
found the N-terminus of vertebrate CavP2a subunit, but notably 
the slowing of inactivation kinetics correlates with the presence of a 
stretch of polybasic residues at the 3' end of the HOOK domain 
common to vertebrate CavPi and CavP^a subunits. While 
alternative splicing of the HOOK domain is primarily responsible 
for promoting the slowing of inactivation of invertebrate CavP 
subunits, alterative splicing of the N-terminus relates to a tissue 
specific expression pattern that is found in CavP subunits. An N- 
terminal "A isoform" composed of exon la/ lb is the predominant 
form expressed in skeletal muscle (although A isoforms of CavP 
subunits are broadly expressed across tissues), while some "B 
isoforms" can have a highly brain-centric expression pattern 
(although this is not true of all "B isoforms" of CavP subunits). 
Specialization such as closed-state inactivation [18,19] are features 
that evolved within vertebrates tailored to specific CavP isoforms 
and Cav2 calcium channel types. Further studies in invertebrates 
will further clarify the nature of the relationship between calcium 
channels and CavP, and it is a model where the numbers of 
possibUities are limited, given that there are only single genes 
coding for Cavl, Cav2, Cav3 and CavP ^ the snail genome. 

Materials and Methods 

Cloning of novel LCavP subunit isoforms 

The calcium channel beta subunit originally cloned and 
expressed from the pond snail Lymnaea stagnalis was LCavPA+ 
and previously described [20], and deposited as GenBank 
Accession # AF484087. Novel LCavPA", LCavPB+ and 
LCavPfi" splice isoforms have been deposited as HM187674. 1, 
HM187675.1 and HM187676.1 respectively. Novel exons in die 
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Figure 7. Snail beta subunit (LCavPA+) splice isoform slow the inactivation kinetics, but has no effect on any other biophysical 
property on snail T-type, LCa„3 channels. (A) Current-voltage relationships (B) Activation and (C) Steady-state inactivation curves. (D) Current 
densities (pA/pF) shown as box plot (box = s.e.m., whisker range = minimum/maximum values). (E) Normalized peak barium current traces illustrated 
as mean, s.e.m. for LCa^S channels co-expressed with LCavP splice isoforms or no Cavp. (F) Curve fitting of inactivation (Tau values) over the steps 
from -110 mV to the range of -70 mV to 0 mV). 
doi:1 0.1 371 /journal.pone.0092941 .g007 



N-terminus and exon 7 were first identified by PGR amplification 
of XZAP cDNA libraries and freshly-isolated genomic DNA. The 
5' end of the LCa^Pa cDNA transcript was extended using 5' 
RACE. Identified sequences were later confirmed by mRNA 
sequences from a published brain transcriptome of Lymnaea 
stagnalis [70], and unpublished genomic sequence spanning the 
LCa^,(3 subunit made available by Daniel Jackson and Angus 
Davison (University of Nottingham). LCa^fi subunit isoforms were 



cloned into mammalian expression vector pMT-2SX(R) using 
PGR inserts with flanking NotI and Xhol restriction sites spanning 
the start and stop codons of the LCa^P gene. DNA primers to 
create the 5' end of the inserts using PGR were designed with a 
standard mammalian KOZ AK sequence downstream of the Not I 
site, GCGGCCGCGAGGATGG, where the underlined ATG 
sequence is the expected start codon. Final expression constructs in 
pMT-2SX(R) vector were confirmed for their expected DNA 
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Figure 8. Snail €3,2 cKiannels nor snail LCBvP subunits do not promote the closed-state inactivation observed for mammalian 
Cav2.2 and CavPs or Ca^Pib subunits. The size of test barium currents relative to the prepulse current after time delays of 0.5, 4, 8, 12, 16, 20 and 
40 ms with coexpression of mammalian Ca^Ps, Ca^pib, CavP2a and snail LCavPA+ with mammalian Ca„2.2 or snail LCav2 calcium channels. A closed- 
state inactivation exhibited where there is an increasing inactivation with increasing time delay, is found only with particular combination of subunits, 
which includes mammalian Ca„2.2 and Ca^Ps or Ca„p,b [18,19]. 
doi:1 0.1 371 /journal.pone.0092941 .g008 



sequences by dideoxy terminator sequencing in both orientations 
(TCAG DNA sequencing facility, The Hospital for Sick Children, 
Toronto, Ontario). 

Sequence comparisons of calcium channel subunits 

Multiple alignments and gene trees of sequences were generated 
in Phylogeny.fr [71]. The running average of similarity across 
Ca^P subunit sequences were generated with PLOTCON in 
EMBOSS [72]. Ca^,|3 subunit sequences from differing exon la/ 
lb, exon 2 and exon 7 splice isoforms in different organisms were 
identified by GenBank BLAST searches. The following are mosdy 
full length Ca^fi subunit sequences used for the generation of 
multiple alignments and gene trees shown in Figure 1 and SI 
{Genus species, GenBank Accession ^): Coanoflagellate [Momsiga 
brevicoUis, XM_00 1748200.1), Poriferan [Amphimedon queenslandica, 
ACUQP 100 1402), Placozoan (Trichoplax adhaerens, ABGPOlO 
00238), Cnidarian [MmatosUlla vectensis, ABAVO 1020487), Lopho- 
trochozoan [Ljmnaea stagnalis, AF484087) and Ecdyzozyan [Dro- 
sophila melanogaster, U11074). The GenBank-derived Ca^l channel 
sequences for Figure ID were: Poriferan {Amphimedon queenslandica, 
XM_003382988), Placozoan {Trichoplax adhaerens, XM_002108 
894), Cnidarian {Mematostella vectensis, XM_00 1639004), Lopho- 
trochozoan {Ljmnaea stagnalis, AF484079) and human Ca„1.2 
(BC 146846). The GenBank-derived Ca^2 channel sequences for 
Figure ID were: Placozoan {Trichoplax adhaerens, XM_002109739), 
Cnidarian {Hydra magnipapillata, XM_0042 10303), Lophotro- 
chozoan {Lymnaea stagnalis, AF484082) and human Ca^2.1 
(AB035727). Genomic and cDNA sequences for human Ca„P 
subunits were derived from NCBI Gene for CACNBl, CACNB2, 
CACNB3, CACNB4 genes. Genomic sequences for snail beta 
subunit were provided by Daniel Jackson and Angus Davison 
(University of Nottingham) and gaps in intron lengths were 



estimated from the genome sequence derived from the closely- 
related air-breathing, freshwater snail Biomphalaria glabrata. 

Calcium channel subunit expression in HEK-293T cells 

Snail LCa^l [15,54] and LCa„2 channels [16,17] from pond 
snail Ljmnaea stagnalis were previously cloned and characterized in 
HEK-293T cells, expressed in pIRES2-EGFP bicystronic vector. 
Mammalian Ca^2.2 (otis) channels were expressed in pMT2 
expression vector, and co-transfected with pTRACER (EGFP). 
Positively-transfected channels were identified by green fluores- 
cence with an Axovert 40 inverted CFL microscope using a GFP 
filter and mercury lamp excitation. Calcium channels were always 
co-transfected with mammalian a^i^ subunit in pMT2 vector, 
with either no Ca^fi subunits, mammalian Ca^fi subunits (CavPn, 
or CavP2a or CavP3) or snaU Ca„P subvmit isoforms. Mammalian 
Ca„P and KjSi subunits were gifts from Gerald Zamponi 
(University of Calgary) and Terry Snutch (University of British 
Columbia). 

Quantitative RT-PCR of calcium channel subunits in snail 
tissues 

mRNA expression was measured using quantitative Real-Time 
PGR (as previously published [14,19,73,74]) using mRNA isolated 
from reproductively-active adult snails from Ljmnaea stagnalis (shell 
lengths of 2.0 to 2.5 cm) illustrated in Figure 4. Adult versus 
sexually immature juvenile mRNA expression was compared in 
Figure 4D, where juveniles have shell lengths of 1 .0 to 1 .5 cm. 
qPCR primers sets (Table SI) were designed to selectively amplily 
LCa^P, LCa^l and LCav2 subunits using a universal primer set, 
and additional primer sets designed to amplify specific exons, such 
as N-terminal exons: LCa^PA and LCa^PB and differing exon 7 
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isoforms: LCovP- and LCavP+. PGR primer specificity was 
confirmed by appropriate sized PGR products amplified fi-om 
pooled template of cDNAs generated fi-om fi-eshly isolated mRNA 
and also cloned cDNAs. PGR primer efficiencies was then 
determined in relative standard curves using 1:5 serial dilutions 
of pooled cDNA (1:5, 1:25, 1:125; and 1:625) as template for real 
time RT-PGR amplification. Triplicate reactions were carried out 
in 96-well PGR plates (Bio-Rad) for each dilution, with each well 
containing 0.5 (xL of serially diluted cDNA, 5 (iL of SsoFastTM 
EvaGreen Supermix (Bio-Rad), 0.5 |XL of each 10 |j,M primer 
from a set, and 3 |J,L of water. PGR amplification, fluorescence 
reading, and melt curve analyses were carried out using a Bio-Rad 
GIOOOTM Thermal Cycler equipped with a GFX96TM Real- 
Time System and run by GFX Manager Software (Bio-Rad). All 
cycle threshold values used for analysis were determined relative to 
the average cycle threshold value of the control gene, HPRTl 
(hypoxanthine phosphoribosyl transferase 1). 

Mammalian HEK-293T cell lines 

We have described our optimized methods for the expression of 
ion channels in mammalian HEK-293T cells and their recording 
using whole-cell patch clamp in an online Jo VE video journal [75] . 
Briefly, HEK-293T are cultured in Dubecco's Modified Eagle's 
Medium (DMEM, Sigma, #D5796) supplemented with 10% Fetal 
Bovine Serum (FBS; Sigma, ^Tli)5l) that had been heat- 
inactivated at 56°G for 30 minutes and 2.5 mL of Penicillin- 
Streptomycin solution (5000 units of penicillin and 5 mg 
streptomycin/mL; Sigma, #P4458). The complete HEK-293T 
media (500 mL) was further supplemented with 5 mL of 100 mM 
sodium pyruvate (Sigma, #88636). In all instances complete 
culture media was heated to 37°G in a water bath before use. 

Subculturing of HEK-293T cells 

To subculture cells, culturing media was removed from the flask 
and cells were washed twice with phosphate buffered saline (PBS; 
6.70 mM KGl, 3.67 mM KH2P04, 10.82 mM Na2HP04- 
7H20, 342 mM NaCl) that had been pre-warmed to 37°G. To 
detach cells, 1 mL of 0.25% Trypsin-EDTA (Sigma, #T4049) 
that had been heated to 37°G was added to the flask. The flask was 
then incubated at 37°G until the majority of cells had detached (to 
a maximum of five minutes). In the meantime, two new flasks were 
fiUed with 5.5 mL of fresh complete media (these will be used for 
the next passage). Also, flasks of cells for transfection were filled 
with 5 mL of fresh complete media. Once the cells were detached, 
5 mL of fresh complete culture media was added and the cells 
were resuspended by pipetting up and down several times. The 
resuspended cells were then divided among the flasks so that each 
flask had a final volume of 6 mL, therefore the flasks to be used in 
the next passage were split 1:12, while the flasks to be transfected 
were spUt 1:6. These cells were then incubated under standard 
conditions. Gells were allowed to grow to 40-50% confluency 
prior to transfection. 

Transient transfection of HEK-293T cells using calcium 
phosphate precipitation 

Transfection of HEK 293T cells was done using a calcium 
phosphate transfection protocol that was carried out by diluting 
4 |i,g of each plasmid to be transfected in 30 |iL of 2.5M GaGl2 
and sterile miUi-Q, water to 300 ^iL. This was then added dropwise 
to 300 hL of 2x HES buffer (280 mM NaGl, 50 mM Hepes, 
1.5 mM Na2HP04-7H20, pH 7.0). The mixture was then mixed 
well and allowed to incubate at room temperature for 20 minutes 
to allow for the formation of calcium phosphate crystals. During 



the incubation, media was removed from the flask to be 
transfected and replaced with 5.4 mL of fresh complete medium. 
After incubation, the calcium phosphate solution was added 
dropwise into the flask of cells. The flask was then incubated under 
standard conditions for eight to 16 hours. After eight to 16 hours, 
the cells were washed twice with PBS pre-heated to 37°C and 
6 mL of fresh complete medium was added. These cells were then 
incubated at 28°G in a humidified 5% G02 atmosphere for two 
days before plating onto glass coverslips. 

Poly-L-Lyslne coating of coverslips and plating cells onto 
coverslips 

Prior to plating HEK cells, sterile, round glass coverslips (Fisher 
Scientific, 12-545-80) were coated in poly-L-lysine. The 
coverslips were spread out in a single layer at the bottom of a 
large (100 mm diameter) culture dish. A dilute poly-L-lysine 
solution was made by diluting 1.5 mL of 0.1% (w/v) poly-L-lysine 
(Sigma, #P8920) in 13.5 mL of miUi-Q^ water. This solution was 
poured over the coverslips and they were allowed to incubate at 
room temperature for one hour. After one hour, the poly-L-lysine 
solution was removed and coverslips were washed twice in 1 5 mL 
of sterile water and then dried in a 56^0 oven for two hours. After 
coating, coverslips were stored at 4°C for up to two weeks. 

Two days after transfection, HEK cells were plated onto glass 
coverslips using the same method. Media was removed from the 
flask and the transfected ( (-Us were washed twice with PBS before 
the addition of 1 mL of tr)'psin, as in section 2.6.3. The cells were 
resuspended in 6 mL of complete media and added to 60 mm 
culture dishes containing the appropriate amount of complete 
media and four to six coverslips so as to give a split ratio of 1 :3 and 
1:4 for cells to be used in patch clamp experiments (HEK cells) 
and 1:6 for cells to be used for antibody staining. These dishes 
were then incubated at standard conditions for three to four hours 
before they were moved to 28°G in a humidified 5% G02 
atmosphere. 

Whole cell patch clamp recording 

Gells were recorded by whole-cell patch clamp method three to 
seven days after plating onto coverslips using an AxoPatch 200B 
amplifier, combined with a Digidata 1440A Data Acquisition 
System and pGLAMP 10 Software (Molecular Devices). Whole 
cell patch clamp recordings were carried out with an Patch 
pipettes for recording with pipette resistances of 2-5 Mfl, and with 
typical access resistance maintained after breakthrough between 4 
and 6 MQ. Only recordings with minimal leak (<\0"/« of peak) 
and small current sizes (<500 pA) in HEK-293T cells were used 
due to loss of voltage clamp above 500 pA. Series resistance was 
compensated to 70% (prediction and correction; lO-^s time lag). 
For all recordings, leak subtraction was preformed offline and data 
was filtered using a 500 Hz Gaussian filter using Glampfit 10.2 
software (Molecular Devices) before further analysis. 

Galcium channel currents were measured in a 20 mM Ba^"*" 
external bath solution: 20 mM BaGlg, 40 mM tetramethylammo- 
nium chloride, 10 mM glucose, 64 mM GsGl, pH 7.2) and patch 
pipettes were filled with internal solution (108 mM cesium 
methane sulfonate, 4 mM MgGlj, 9 mM HEPES, 9 mM ethylene 
glycol tetraacetic acid, pH 7.2). During all recording sessions, cells 
were maintained at a holding potential of — 1 00 mV. Prior to 
recording, a test step to peak current (10 mV for LGa^l or 30 mV 
for LGav2) was performed to indicate if the cell was producing a 
current and if the state of the patch was suitable for gathering data. 
All data derived from patch clamp recording was statistically 
analyzed using one-way analysis of variation (ANOVA) tests online 
at http://www.danielsoper.com/statcalc3. 
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Electrophysiology protocols for generating current- 
voltage relationships and steady-state inactivation 

The current-voltage (IV) relationship was assessed by stepping 
cells from — 100 mV to —50 mV for 450 ms and then increasing 
the voltage step by 10 mV in successive sweeps until reaching a 
potential of 50 mV. To analyze the effect of LCavP subunits on the 
steady-state inactivation of the channels cells were stepped from — 
100 mV to peak voltage for 150 ms, then allowed to recover for 
1 s before being held at — 1 00 mV untU all channels had reached 
inactivation (this conditioning step may take up to 15 seconds) 
then immediately stepping to peak again for 150 ms. In successive 
sweeps, the conditioning potential used to inactivate channels was 
increased by 10 mV each time untU reaching 30 mV. 

Analyses of activation and current-voltage curves 

To create activation and IV plots, current (I in pA) at each 
voltage step was normalized to peak current (I„,ax in pA) for that 
cell. Each normalized IV relationship was plotted and the 
individual reversal potentials (E^ev; in mV) were determined by 
calculating the y-intercept of the linear portion of the IV curve (+ 
20 to -1-40 for LCa^l ; -1-30 to -1-50 for LCa^,2). The reversal potential 
represents the membrane potential at which the driving force for 
Ca^"*" (or Ba~^) influx is equal to the driving force for Ca^"*" (or Ba^"*^ 
efflux and there is no net movement of Ca^^ (or Ba^""^ through the 
calcium channels. The conductance (rate of ions flowing through 
the channel) was then determined for each voltage step using the 
equation: G = (I/Imax)/(V— E^ev), where G represents conductance 
(pS) of the channel and V represents the test voltage (mV). For 
each trace maximum conductance (G„-j,J was then determined. 
The mean and standard error of the mean (s.e.m.) were calculated 

for (I/Imax), Erev and Gmax- 

Activation plots were created using Origin 8 software 
(OriginLab, Northampton, MA) by importing IV data and then 
performing a Boltzmann transformation using the equation: 
Conductance (gj = l/(V—'E^„). The Boltzmann-transformed data 
was then normalized and then a scatter of normalized activation 
versus voltage was created and curve fitted with the following 
Boltzmann equation: G/Gniax= l/(l-l-e((V— Vo.5)/Ka)), where 
Vo.5 and represent the half-activation voltage and slope factor 
of the activation curves, respectively. The half-activation potential 
is the voltage at which 50% of channels is expected to be open. 
The slope factor of the activation curve is a value that represents 
the delay of channels transitioning from the closed to active (open) 
position. The normalized activation data [(I/Imax)) Vq.s, KJ were 
then averaged and the SEM was determined. The mean data was 
then plotted and a curve was simulated using the fitted parameters. 
IV plots were created by plotting voltage against mean (I/In,ax) — 
SEM and then fitted with an Ohmic-Boltzmann curve using the 
following equation: I/Imax = G„iax(V— Erev)/(l+e^ vo.s^Ka^ 

Analyses of steady-state inactivation 

Steady-state inactivation values were determined by dividing the 
current after the inactivating (conditioning) pulse by the current 
prior to the inactivating pulse. Each data set was then imported 
into Origin and curve fit with the Boltzmann equation: (1/ 
Imiix) = 1/(1+"*^ vo.5)/Ki^ w'here Vo.s and K; represent the half- 
inactivation voltage and the slope factor of inactivation, respec- 
tively. The half-inactivation potential is the voltage at which 50% 
of the channels are available and the other 50% of channels have 
transitioned from the open to the inactivated state. The values of 
V().5 and K; (determined using Origin 8) and (I/Imax) were then 
averaged and the standard error of each means were calculated. 
The mean steady-state activation data was then plotted ± SEM 



and this data was used to simulate a Boltzmann curve derived 
from the fitted parameters. The relative inactivation of LCa,,l and 
LCav2 currents were compared as the fraction of current 
remaining after 250 ms (R250). 

Western blotting of LCavP subunits 

LCa^P protein was identified on Western blots with rabbit anti- 
LCa^P subunit antibody made against KLH-coupled synthetic 
peptide with LCa„P sequence: SLDEEKEALRRET, which is 
downstream of the N-terminus and common to all LCa^P isoforms 
[20]. Western blots were prepared from protein homogenates 
isolated from HEK cell lysates, five days post-transfection. Total 
protein was separated by SDS-PAGE, and the contents of the gels 
were transferred onto 0.45 [>M nitrocellulose membrane (Mandel, 
# W- 1040 1196) overnight at 90 mA at 30 V. Before applying 
antibody, the nitrocellulose membrane must be blocked to reduce 
nonspecific interactions. Blocking buffer [10 mM Tris, 100 mM 
NaCl, 0.1% Tween-20(v/v), 5%(w/v) skim milk powder] was 
added to the membrane in a shallow dish (enough to completely 
cover the nitrocellulose), placed on a shaker and allowed to 
incubate at room temperature for two hours. 

After blocking, the membrane was washed five times with 
Tween-Tris-BuflFered Saline [TTBS; 10 mM Tris, 100 mM NaCl, 
0. 1 % Tween-20(v/v)] for five minutes at room temperature. Next, 
the primary antibody (rabbit anti-LCavP) was diluted 1:2000 in 
blocking buffer and poured over the membrane after the last wash 
with TTBS. The primary antibody was incubated with the 
membrane overnight at 4°C. The following morning, the antibody 
was removed and the membrane was washed five times in TTBS 
for five minutes at room temperature, and then once in blocking 
buffer for 15 minutes at room temperature. After the second 
block, the secondary antibody (Goat anti-rabbit IgG coupled to 
horseradish peroxidase; Invitrogen, 77^65-6120) was diluted 1:5000 
in blocking buffer and added to the dish containing the 
membrane. The secondary antibody was incubated with the 
membrane for 30 minutes at room temperature before it was 
decanted. The membrane was then washed five times in TTBS for 
five minutes at room temperature. The membrane was then ready 
for chemiluminescent staining to detect the presence of LCavP 
subunits. 

LCa^P subunits were detected using the chemiluminescent stain 
ECL which was made by combining 20 mL of freshly prepared 
solution one (2.5 mM luminol, 396 |i.M p-coumaric acid, 100 mM 
Tris, pH 8.5) and 20 mL of freshly prepared solution two 
(100 mM Tris, 12 uL of 30% H2O2, pH 8.5). Solutions one and 
two were combined, poured over the membrane and incubated for 
one minute at room temperature. The blots were then exposed to 
x-ray ftlm (Kodak, #819 4540) for one minute in a darkroom 
before development using an automated developer. Films were 
then compared to nitrocellulose membranes and the ladder was 
marked onto the film by hand. 

Potential caveats in studies of LCavp subunits expressed 
in HEK-293T cells 

We observe that snail LCa^P subunit isoforms promote 
membrane expression of calcium channels in HEK-293T cells, 
but the biophysical effects are relatively minor compared to the 
mammalian CavP subunits. A potential concern is a weak signal to 
background ratio perhaps due to a lack of saturation of 
heterologously-expressed Ca^P subunits in every cell recording 
that is also to contain calcium channels and a co-expressed agS 
subunit in a standard calcium-phosphate transfection protocol. 
Another issue is potentially contaminating, low levels of human 
calcium channels subunits reported in HEK-293T cells [76]. A 
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fetal form of CavPS subunit is detectable in Xenopus oocytes 
[77,78], but little to no endogenous CavP subunits have been 
reported in HEK-293T cells [79] . We are confident that our data 
represents the snail Ca^fi subunit isoforms effects on expressed 
snail calcium channels. Positively-expressing LCa^l and LCa^2 
HEK-293T cells corresponded largely to the green EGFP 
fluorescence of transfected cells, since EGFP and the calcium 
channels were derived from the same mRNA on bicistronic vector 
pIRES2-EGFP. LCa^fiA and LCa^Pa subunit expression can be 
identified by anti-rabbit LCa^P specific antibody in Western blots 
of transfec ted HEK cell lysates at appropriate protein sizes 
(Figure S3). 

Expression of mammalian Qtj-S subunits can cause dramatic 
shifts in voltage-sensitivities [7,80], and perhaps the full contribu- 
tion of moUuscan Ca^P subunits will not be completely described 
without co-expression of native a2-S subunits. Mining of moUuscan 
genomes reveals three <x.2-b subunits: a generalized one that 
resembles a2-5i and a2-52 subunit in structure, a second brain 
specific one that is a homolog to Drosophila straight] acket [81] and 
human ocg'^s subunits, and a third more invertebrate specific 0(2-5 
subunit. 

Supporting Information 

Figure SI Multiple alignment of Cavp subunits includ- 
ing those from a cnidarian (Nematostella), nematode 
(Caenorhabditis), ecdysozoan (Drosophila), lophotro- 
chozoan (Lymnaea) and human gene isoforms (Gavpib, 
CavP2a, CavpS, CavP4b). Highly conserved SH3 and GK 
domains are illustrated, also conserved secondary structures (a 
hehces and P sheets), and calcium channel (AID) binding residues 
(blue residues) reported in crystal structures of CavP subunits. 
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Exon boundaries (red lines) are indicated. Red boxes surrounding 
a base indicates that the intron splits between an amino acid base. 
(TIF) 

Figure S2 Genomic sequences spanning exon 7 illus- 
trating the alternate acceptor sites in Lymnaea snail 
calcium channel beta subunit, (LCavpf and LCavp— ) 
which generates a seven amino acid optional exon (7A-I-) 
or not exon (7A— ). Note that skipping of exon 7 generates a 
change in reading frame and truncated LCavP. The frame shift 
occurs because exon 6 ends within a codon after the first 
nucleotide (phase 1), whereas the intron preceding exon 8 is 
located between codons (phase 0). 
(TIF) 

Figure S3 Snail LCavP - specific rabbit antibody local- 
izes LCavpA+ and LCavPB+ proteins in transfected HEK 
ceUs by Western blotting. SnaU LCa,2 in pIRES2-EGFP 

vector were transfected alone (mock) or with coexpressed 
LCavPA-l- and LCavPg-H plasmids in HEK cells. Expressed 
LCavPA-l- and LCavPs-l- proteins are identifiable on Western blots 
of transfected HEK cell lysates at appropriate size (62.8 kDa and 
62.6 kDa), respectively. 
(TIF) 

Table SI Q,uantitative Real Time PCR (qPCR) primer 
sequence parameters used in Figure 4. 

(DOCX) 

Author Contributions 

Conceived and designed the experiments: TFD ANB AS JDS. Performed 
the experiments: TFD ANB AS JP ST JDS. Analyzed the data: TFD ANB 
AS JP ST DJ .AD JDS. Contributed reagents/materials/analysis tools: DJ 
AD JDS. Wrote the paper: TFD ANB AS JDS. 



17. Huang X, Scnatorr A, Dawson IF, Quan Spaliord JD (2010) G-proteins 
modulate invertebrate synaptic calcium channel (LCav2) differently from the 
classicEd voltage-dependent regulation of mammalian Cav2.1 and Cav2.2 
channels. J Exp Biol 213: 2094-2103. 

18. Senatore A, Spafford JD (2010) Transient and big are key features of an 
invertebrate T-type channel (LCav3) from the central nervous system of 
Lymnaea stagnalis. J Biol Ghcm 285: 7447-7458. 

19. Senatore A, Spafibrd JD (2012) Gene transcription and spHcing of T-type 
channels arc cvoliilionarily-conserved strategies for regulating channel expres- 
sion and gating. PLoS One 7: e37409. 

20. Spafford JD, van MJ, Larsen P, Smit AB, Syed Nl, et al. (2004) Uncoupling of 
calcium channel alphal and beta subunits in developing neurons. J Biol Chem 
279: 41157-^1167. 

21. Salvador-Recatala V, Schneider T, Greenberg RM (2008) Atypical properties of 
a conventional calcium channel beta subunit from the platyhelminth 
Schistosoma mansoni. BMC Physiol 8: 6. 

22. Kimura T, Kubo T (2003) Cloning and functional characterization of squid 
voltage-dependent Ca2+ channel beta subunits: involvement of N-terminal 
sequences in differential modulation of the current. Neurosci Res 46: 105-1 17. 

23. Gens T, Rousset M, GoUet G, Raymond V, Dcmares f, et al. (2013) 
Gharacterization of the first honeybee Ca(2)(+) channel subunit reveals two 
novel species- and splicing-speeific modes of regulation of channel inactivation. 
Pllugcrs Arch 465: 985-996. 

24. Chien AJ, Carr KM, Shirokov RE, Rios E, Hosey MM (1996 Idcniilit ation of 
palmitoylation sites within the L-type calcium channel beta2a subunit and effects 
on channel function. J Biol Chem 271: 26465-26468. 

25. Qin N, Platano D, Olcese R, Costantin JL, Stefani E, et al. (1998) Unique 
regulatory properties of the type 2a Ca2+ channel beta subunit caused by 
palmitoylation. Proc Natl Acad Sci USA 95: 4690-4695. 

26. Restituito S, Cens T, Barrere C, Geib S, Galas S, et al. (2000) The [beta] 2a 
subunit is a molecular groom for the Ca2+ channel inactivation gate. J Neurosci 
20: 9046-9052. 

27. Yasuda T, Lewis RJ, Adams DJ (2004) Overexprcs.sed C.:a!\ ;beia3 inhibits N- 
tyjie (Gav2.2) calcium channel currents through a hyperpolarizing shift of ultra- 
slow and closed-state inactivation. J Gen Physiol 123: 401-416. 

28. Patil PG, Brody DL, Yue DT (1998) Preferential closed-state inactivation of 
neuronal calcium channels. Neuron 20: 1027—1038. 



PLOS ONE I www.plosone.org 



16 



April 2014 | Volume 9 | Issue 4 | e92941 



Invertebrate Cavp Regulation of Calcium Channels 



29. Hanlon MR, Bcrrow NS, Dolphin AC, Wallace BA (1999) Modelling of a 
voltage-dependent Ca2+ channel beta subunit as a basis for understanding its 
functional properties. FEBS Lett 445: 366-370. 

30. te Velthuis AJ, Admiraal JF, Bagowski CP (2007) Molecular evolution of the 
MAGUK family in metazoan genomes. BMC Evol Biol 7: 129. 

31. Chen YH, Li MH, Zhang Y, He LL, Yamada Y, et al. (2004) Structural basis of 
the alphal-beta subunit interaction of voltage-gated Ca2-I- channels. Nature 429: 
675-680. 

32. Opatowsky Y, Chen CC ('ampbell KP, Hirsch JA (2004) Structural analysis of 
the voltage-dependent calcium channel lieta siiliiinil liaiclional core and its 
complex with the alpha 1 interaction domain. iVcuron 42: 387-399. 

33. Van PF, Clark K.\, Chatelain FC, Minor DLJr (2004) Structure of a complex 
between a voltage-gated calcium channel beta-subunit and an alpha-subunit 
domain. Nature 429: 671-675. 

34. De WM, Witcher DR, PragncU M, Liu H, Campbell KP (1995) Properties of the 
alpha 1-beta anchoring site in voltage-dependent Ca2+ channels. J Biol Chem 
270: 12056-12064. 

35. Ebert AM, McAnelly CA. Handschy AV, Mueller RL, Home WA, et al. (2008) 
(Genomic organization, expression, and phylogenctic analysis of Ga2-I- channel 
beta4 genes in 13 vertebrate species. Physiol (ienomics 35: 133—144. 

36. Ebert AM, McAnelly CA, Srinivasan A, Mueller RL, Garrity DB, et al. (2008) 
The calcium channel beta2 (CACNB2) subunit repertoire in teleosts. BMC Mol 
Biol 9: 38. 

37. Shepard S, McCrcary M, Fcdorov A (2009) The peculiarities of large intron 
splicing in animals. PLoS One 4: e7853. 

38. Currie KP (2010) G protein modulation of CaV2 voltage-gated calcium 
channels. Channels (Austin) 4: 497-509. 

39. Obermair GJ, Tuluc P, Flucher BE (2008) Auxiliary Ca(2+) channel subunits: 
lessons learned from muscle. Curr Opin Pharmacol 8: 311—318. 

40. Schredelseker J, Di BV, Obermair GJ, Felder ET, Flucher BE, et al. (2005) The 
beta la subunit is essential for the assembly of dihydropyridine-receptor arrays in 
skeletal muscle. Proc Natl Acad Sci USA 102: 17219-17224. 

41. Murakami M, Nakagawasai O, Yanai K, Nunoki K, Tan-No K, et al. (2007) 
Modified behavioral characteristics following ablation of the voltage-dependent 
calcium channel beta3 subunit. Brain Res 1160: 102—112. 

42. Witcher DR, De WM, Sakamoto J, Franzini-Armstrong C, Pragnell M, et al. 
(1993) Subunit identification and reconstitution of the N-type Ca2+ channel 
complex purified from brain. Science 261: 486-489. 

43. Miranda-Laferte E, Schmidt S, Jara AC, Neely A, Hidalgo P (2012) A short 
polybasic segment between the two conserved domains of the beta2a-subumt 
modulates the rate of inactivation of R-type calcium channel. J Bio! Chem 287: 
32588-32597. 

44. Cohen RM, Foell JD, Balijepalli RC, Shah V, Hell JW, et al. (2005) Unique 
modulation of L-type Ca2-I- channels by short auxiliary beta Id subunit present in 
cardiac muscle. AmJ Physiol Heart Circ Physiol 288: H2363-H2374. 

45. Foell JD, Balijepalli RC, Delisle BP, Yunker AM, Robia SL, et al. (2004) 
Molecular heterogeneity of calcium channel beta-subunits in canine and human 
heart: evidence for differential subcellular localization. Physiol Genomics 17: 
183-200. 

46. Harry JB, Kobrinsky E, Abemethy DR, Soldatov NM (2004) New short splice 
variants of the human cardiac Cavbeta2 subunit: redefining the major functional 
motifs implemented in modulation of the Cavl.2 channel. J Biol Chem 279: 

46367-46372. 

47. Hibino H, Pironkova R, Onwumerc O, Roussct M, Charnct P, ct al. (2003) 
Direct interaction with a nuclear protein and regulation of gene silencing by a 
variant of the Ca2-l-channel beta 4 subunit. Proc Natl Acad Sci USA 100: 307- 
312. 

48. Subramanyam P, Obermair GJ, Baumgartner S, Gebhart M, Striessn% J, et al. 
(2009) Activity and calcium regulate nuclear targeting of the calcium channel 
beta4b subunit in nerve and muscle cells. Channels (Austin) 3: 343-355. 

49. Xu X, Lee YJ, Holm JB, Terry MD, Oswald RE, et al. (201 1) The Ca2+ 
channel beta4c subunit interacts with heteroehromatin protein 1 via a PXVXL 
binding motif. J Biol Chem 286: 9677 9687. 

50. Zhang Y, Yamada Y, Fan M, Bangaru SD, Lm B, et al. (2010) The beta subunit 
of voltage-gated Ca2-I- channels interacts with and regulates the activity of a 
novel isoform of Pax6. J Biol Chem 285: 2527-2536. 

51. Jangsangthong W, Kuzmenkina E, Khan IF, Matthes J, Hullin R, et al. (2010) 
Inactivation of L-type calcium channels is determined by the length of the N 
terminus of mutant beta(l) subunits. Pflugers Arch 459: 399-411. 

52. Herzig S, Klian IF, Grundemann D, Matthes J, Ludwig A, et al. (2007) 
Mechanism of Ca(v)1.2 channel modulation by the amino terminus of cardiac 
beta2-subunits. FASEB J 21: 1527-1538. 

53. Richards MW, Leroy J, Pratt WS, Dolphin AC (2007) The HOOK-domain 
between the SH3 and the GK domains of Cavbeta subunits contains key 
determinants controUing calcium channel inactivation. Channels (Austin) 1 : 92— 
101. 



54. Spafford JD, Dunn T, Smit AB, Syed NI, Zamponi GW (2006) In vitro 
characterization of L-type calcium channels and their contribution to firing 
behavior in invertebrate respiratory neurons. J Neurophysiol 95: 42-52- 

55. Christel C, Lee A (2012) Ca2-Kdependent modulation of voltage-gated Ca2-I- 
channels. Biochim Biophys Acta 1820: 1243-1252. 

56. Almagor L, Chomsky-Hecht O, Ben-Mocha A, Hendin-Barak D, Dascal N, 
et al. (2012) The role of a voltage-dependent Ca2-I- channel intracellular linker: a 
structure-flinction analysis. J N('urosci 32: 7602—7613. 

57. Vitko I, Shcheglovito\' A. Baumgart }P, Arias-Olguin II, Murbartian J, et al. 
(2008) Orientation of the (^aleium ('hannel beta relative to the alpha(l)2.2 subunit 
is critical for its regulation of channel activity. PLoS One 3: e3560. 

58. Zhang Y, Chen YH, Bangaru SD, He L, Abele K, et al. (2008) Origin of the 
voltage dependence of G-protein regulation of P/Qjtype Ca2+ channels. 
J Neurosci 28: 14176-14188. 

59. Bae J, Suh EJ, Lee C (2010) Interaction of T-type calcium channel Ca(V)3.3 
with the beta-subunit. Mol Cells 30: 185-191. 

60. Duliel SJ, Altier G, Chaumont S, Lor\'- P, Bourinet E, et al. (2004) Plasma 
membrane expression of T-t)pe calcium channel alpha(l) subunits is modulated 
by high voltage -activated auxiliary subunits. J Biol Chem 279: 29263—29269. 

61. Laccrda AE, Perez-Reyes E, Wei X, Castellano A, Brown AM (1994) T-type 
and N-type calcium channels of Xenopus oocytes: evidence for specific 
interactions with beta subunits. Biophys J 66: 1833-1843. 

62. Arias JM, Murbartian J, Vitko I, Lee JH, Perez-Reyes E (2005) Transfer of beta 
subunit regulation from high to low voltage-gated Ca2-I- channels. FEBS Lett 
579: 3907-3912. 

63. Lambert RC, Maulet Y, MoutonJ, Beattie R, Volsen S, et al. (1997) T-t>pe 
Ca2-I- current properties are not modified by Ga2-I- channel beta subunit 
depletion in nodosus ganglion neurons. J Neurosci 17: 6621—6628. 

64. Leurangucr V, Bourinet E, Lory P, NargeotJ (1998) Antisense depletion of beta- 
subunits fails to affect T-type calcium channels properties in a neuroblastoma 
cell Hne. Neuropharmacology 37: 701—708. 

65. Senatore A, Zhorov BS, Spafford JD (2012) Gav3 T-type calcium channels. 
WIREs Membr Transp Signal 1: 467-491. 

66. Fang K, Colecraft HM (2011) Mechanism of auxiliary beta-subunit-mediated 
membrane targeting of L-t\'jie (Ca(V)1.2) channels. J Physiol 589: 4437-4455. 

67. Altier C, Garcia-Caballero A, Simms B, You H, Chen L, et al. (2011) The 
Cavbeta subunit prevents RFP2-mediated ubiquitination and proteasomal 
degradation of L-type channels. Nat Neurosci 14: 173-180. 

68. Rougier JS, Albesa M, Abriel H, Viard P (2011) Neuronal precursor cell- 
expressed developmentaUy down-regulated 4-1 (NEDD4-1) controls the sorting 
of newly synthesized Ca(V)1.2 calcium channels. J Biol Chem 286: 8829-8838. 

69. Waithe D, Ferron L, Page KM, Ghaggar K, Dolphin AC (201 1) Beta-subunits 
promote the expression of Ca(\')2.2 channels by reducing their proteasomal 
degradation. J Biol Chem 286: 9598-961 1. 

70. Sadamoto H, Takahashi H, Okada T, Kenmoku H, I'oyota M, et al. (2012) De 
novo sequencing and transcriptome analysis of the central nervous system of 
mollusc Lymnaea stagnalis by deep RNA sequencing. PLoS One 7: e42546. 

71. Dereeper A, Guignon V, Blanc G, Audic S, Buffet S, et al. (2008) Phylogeny.ft: 
robust phylogenetic analysis for the non-specialist. Nucleic Acids Res 36: W465— 
W469. 

72. Rice P, Longden I, Bleasby A (2000) EMBOSS: the European Molecular 
Biology Open Software Suite. Trends Genet 16: 276—277. 

73. Senatore A, Spafford JD (2013) A uniquely adaptable pore is consistent with 
NALCN being an ion sensor. Chaimcls (Austin) 7. 

74. Senatore A, Guan W, Boone AN, Spafford JD (2014) T-type calcium channels 
become sodium chEinnels using alternate extracellular turret residues outside the 
selectivity filter. Journal of Biological Chemistry. 

75. Senatore A, Boone AN, Spafford JD (2011) Optimized transfection strategy for 
expression and electrophysiological recording of recombinant voltage-gated ion 
channels m HEK-293T cells. J Vis Exp. 

76. Berjukow S, Doring F, Froschmayr M, Grabner M, Glossmann H, et al. (1996) 
Endogenous calcium channels in human embryonic kidney (HEK293) cells. 
Br J Pharmacol 118: 748-754. 

77. Canti C, Davies A, Bcrrow NS, Butcher AJ, Page KM, et al. (2001) Evidence for 
two concentration-dependent processes for beta-subunit effects on alphalB 
calcium channels. BiophysJ 81: 1439-1451. 

78. Tareilus E, Roux M, Qin N, Olcese R, Zhou J, et al. (1997) A Xenopus oocyte 
beta subunit: evidence for a role in the assembly/ expression of voltage-gated 
calcium channels that is separate from its role as a regulatory subunit. Proc Nad 
Acad Sci USA 94: 1703-1708. 

79. Yasuda T, Chen L, Barr W, McRory JE, Lewis RJ, et al. (2004) Auxiliary 
subunit regulation of high-\'oltage activated calcium channels expressed in 
mammalian cells. Eur J Neurosci 20: 1—13. 

80. Dolphin AC (2012) The alpha(2)delta subunits of voltage-gated calcium 
channels. Biochim Biophys Acta. 

81. Ly CV, Yao CK, Vcrstrckcn P, Ohyama T, Bellcn HJ (2008) straightjacket is 
required for the synaptic stabilization of cacophony, a voltage-gated calcium 
channel alphal subunit. J CeU Biol 181: 157-170. 



PLOS ONE I www.plosone.org 



17 



April 2014 | Volume 9 | Issue 4 | e92941 



